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CN ■ ABSTRACT 

\o ■ 

, Aims. The large majority of BL Lacertae objects belonging to the 1 Jy sample, the class prototype for radio-selected sources, are thought to 
• emit most of their synchrotron power in the far infrared band. Ironically, this spectral region is very sparsely sampled, with only a minority of 
' the objects having IRAS data (most of them being upper limits or low-quality detections). We aim at filling this infrared gap by presenting new, 
simultaneous ISOCAM and ISOPHOT observations over the 7 - 200yum range (observer's frame) for half the sample. A precise measurement 

^ i' of the position of the synchrotron peak frequency, Vpeak, can provide direct information about particle acceleration mechanisms and constrain 

JL the inverse Compton radiation that will be detected by up-coming new y-ray missions. 

^ Methods. We have observed seventeen 1 Jy BL Lacertae objects with the camera and the photometer on board the Infrared Space Observatory 
(ISO) satellite. Given the intrinsic variability of these sources, the data were taken by concatenating the pointings to ensure simultaneity. The 
ISOPHOT data reduction was done employing a novel correction, which mitigates the effect of chopping for faint sources. 

IL" Results. Using our new ISO data, complemented by nearly-simultaneous radio and optical observations for ten and four objects respectively, 
. 5^ and other multi-frequency data, we have built the spectral energy distributions of our sources (plus a previously published one) and derived 
the rest-frame Vpeaj^. Its distribution is centered at ~ 10'^ Hz (~ 30yum) and is very narrow, with ~ 60% of the BL Lacs in the 1 - 3 x 10'^ 

5_j Hz range. Given our set of simultaneous infrared data, these represent the best determinations available of the synchrotron peak frequencies 
for low-energy peaked BL Lacs. A comparison with previous such estimates, based on non-simultaneous optical and near infrared data, may 
indicate strong Vp^^k variations in a number of sources, possibly associated with large flares as observed in the high-energy peaked BL Lac 
MKN 501. 



Key words. BL Lacertae objects ~ Infrared: galaxies - Methods: data analysis - Radiation mechanisms: non-thermal 



1. Introduction 

BL Lacertae objects constitute one of the most extreme classes 
of active galactic nuclei (AGN), distinguished by their high lu- 
minosity, rapid variability, high (> 3%) optical polarization, ra- 
dio core-dominance, apparent superluminal speeds, and almost 
complete lack of emission lines (e.g., UiTy & Padovani 1995 1. 
The broad-band emission in these objects, which extends from 
the radio to the gamma-ray band, appears to be dominated by 
non-thermal processes from the heart of the AGN, undiluted by 
the thermal emission present in other AGN. Therefore, BL Lacs 
represent the ideal class to study to further our understanding 
of non-thermal emission in AGN. 
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Synchrotron emission combined with inverse Compton 
scattering is generally thought to be the mechanism responsible 
for the production of radiation over such a wide energy range 
(e.g., Ghisellini et al. .1998J . The frequency at which most of 
the synchrotron power is emitted, Vpeak, ranges across several 
orders of magnitude, going from the far-infrared to the hard 
X-ray band. Sources at the extremes of this wide distribution 
are referred to as low-energy peaked (LBL) and high-energy 
peaked (HBL) BL Lacs, respectively (Giommi & Padovani 
il994t Padovani & Giommi . 1995j . Radio-selected samples in- 
clude mostly objects of the LBL type, while X-ray selected 
samples are mostly made up of HBL. 

It is ironic that the spectrum of LBL in the infrared (IR) 
band, where these sources' output is thought to be largest, is 
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Table 1. Sample properties. 
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Name 


RA(J2000) 


Dec(J2000) 


z 




FsGHz 












Jy 


PKS 0048-097 


00 50 41.3 


-09 29 06 


>0.2 


16.5 


2.0 


PKS 0118-272 


01 20 31.7 


-27 01 24 


>0.559 


17.0 


1.1 


PKS 0138-097 


01 41 25.8 


-09 28 44 


0.733 


17.0 


1.2 


PKS 0426-380 


04 28 40.4 


-37 56 19 


1.110 


18.0 


1.2 


S5 0454+844 


05 08 42.4 


+84 32 05 


> 1.340 


19.0 


1.4 


PKS 0537-441 


05 38 50.4 


-44 05 09 


0.896 


13.5 


4.0 


PKS 0735+178 


07 38 07.1 


+ 17 42 27 


>0.424 


14.5 


2.0 


PKS 1144-379 


11 47 01.4 


-38 12 10 


1.048 


16.5 


1.6 


OQ530 


14 19 46.6 


+54 23 15 


0.152 


14.5 


1.1 


PKS 1519-273 


15 22 37.7 


-27 30 10 


>0.2 


18.5 


2.4 


PKS 1749+701 


17 48 32.8 


+70 05 51 


0.770 


16.5 


1.1 


PKS 1749+096 


17 51 32.8 


+09 39 01 


0.320 


15.5 


1.9 


S5 1803+784 


18 00 45.7 


+78 28 04 


0.684 


16.5 


2.6 


4C 56.27 


18 24 07.2 


+56 51 00 


0.664 


18.5 


1.7 


PKS 2131-021 


21 34 10.3 


-01 53 17 


1.285 


20.0 


2.1 


PKS 2240-260 


22 43 26.4 


-25 44 31 


0.774 


17.5 


1.0 


PKS 2254+074 


22 57 17.3 


+07 43 12 


0.190 


17.0 


1.2 



very sparsely known (e.g., Sambruna et al.'1996"; Giommi et al. 
[2002. Massaro et al. 2QQ5 l. IR data for the 1 Jy, the prototype 
radio-selected BL Lac sample are, for example, very sparse. 
The purpose of this paper is to remedy this situation by present- 
ing ISOCAM and ISOPHOT observations of 1 Jy BL Lacs cov- 
ering the 6 .7 - 200//m range. Since Vpeak Tpeak^^' where ypeak 
is the Lorentz factor of the electrons emitting most of the radi- 
ation, 6 is the Doppler factor, and B is the strength of the mag- 
netic field, its precise measurement can provide direct infor- 
mation about particle acceleration mechanisms. The position 
of Vpeak is also important for the up-coming Astro-rivelatore 
Gamma a Immagini Leggero (AGILE') and Gamma ray Large 
Area Space Telescope (GLAST^) missions, as they will sample 
the y-ray band close to the peak of the inverse Compton emis- 
sion from LBL and flat-spectrum radio quasars, whose location 
depends on that of the synchrotron one. 

In Sect. 2 we present our sample. Sect. 3 discusses the 
Infrared Space Observatory (ISO) observations, data analysis, 
and results, while in Sect. 4 we derive the spectral energy dis- 
tributions and Vpeak values. Finally, Sect. 5 summarizes our 
conclusions. Throughout this paper spectral indices are written 

5v oc v^" . 

2. The sample 

The 1 Jy sample of BL Lacs is the only sizeable, complete sam- 
ple of radio bright BL Lacs. It includes 34 objects with radio 
flux > 1 Jy at 5 GHz (Stickel et al. TW). AU 1 Jy BL Lacs 
have been studied in detail in the radio and optical bands; all 
objects have also soft X-ray data, primarily from ROSAT (Urry 
et al. \199ES. Most of the sources have also been detected by 
the Wilkinson Microwave Anisotropy Probe (WMAP) satellite 
(Bennett et al. I2003t and about 40% of them have also hard 

' |]http : //agile .rm.iasf rcnr. it/I 
^ pittp : //glast . gsfc . nas^gov/] 



X-ray (BeppoSAX) data (e.g., Padovani et al. 120041 and refer- 
ences therein). In the IR band, on the other hand, 1 Jy BL Lacs 
are not very well studied. For example, only ~ 40% of 1 Jy 
BL Lacs have IRAS data, and many of these are of low qual- 
ity or only upper limits. We selected for ISO observations all 
1 Jy BL Lacs visible by the satellite and not already part of 
other programs. This included seventeen 1 Jy BL Lacs (or 50% 
of the sample). All sources are LBL. The object list and basic 
characteristics are given in Tab.^ Accurate radio positions and 
redshifts come from the NASA/IPAC Extragalactic Database 
(NED), mean R magnitudes are from Heidt & Wagner (,1996.) . 
while 5 GHz radio fluxes are from Stickel et al. J1991> . 

3. Observations, data analysis, and results 

The ISO satellite (Kessler et al. ll996> was equipped with a 60 
cm Ritchey-Chretien telescope and had four instruments. For 
our observations we used the camera (ISOCAM) and the pho- 
tometer (ISOPHOT). Given the intrinsic variability of BL Lacs, 
we needed to observe each object at a single epoch. Thus, the 
three different Astronomical Observation Templates (AOTs) 
we used (see below) were concatenated in order to warrant as 
much simultaneity as possible for the measurements. 

3.1. ISOCAM 

All sources were observed with the 32 x 32 pixel ISOCAM 
(Blommaert et al. l2003> in staring mode (AOT = CAMOl) us- 
ing the long-wavelength (LW) detector in three different bands: 
LW2 (A = 6.1 fim), LW7 (A = 9.6;um), and LW3 (A = 14.3;um). 
The pixel scale was 6" and the total field of view 196 x 196 
arcsec". The total exposure times were computed as Tim x A^exp, 
where Texp = 2.1s is the integration time of a single exposure 
and A^exp is the number of exposures. Typical exposure times 
are 100 s, 100 s, and 75 s for the LW2, LW7, and LW3 filters, re- 
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Table 2. ISO journal of observations. 



Table 3. ISOCAM fluxes" and slopes. 



Name 


ISOJd 


Observing date 


Name 


LW2 


LW7 


LW3 


ay 




CAM/P2/C100&C200 






[6.7;um] 


[9.6;um] 


[14.3yum] 




PKS 0048-097 


38701401/02/03 


1996 Dec 8 


PKS 0048-097 


18.2±2.0 


19.2+6.8 


25.7+5.7 


0.4+0.4 


PKS 0118-272 


74400104/05/06 


1997 Nov 28 


PKS 0118-272 


18.2+2.7 


20.1+7.3 


35.7+5.6 


0.9+0.3 


PKS 0138-097 


76301007/08/09 


1997 Dec 17 


PKS 0138-097 


9.1±2.2 


11.9+5.5 


18.9+5.2 


1.0+0.6 


PKS 0426-380 


69702410/11/12 


1997 Oct 12 


PKS 0426-380 


2.5±1.2 


<14.7 


19.0+5.6 


2.7+0.7 


S5 0454+844 


68700413/14/15 


1997 Oct 2 


S5 0454+844 


5.5±2.2 


<16.2 


13.5+5.7 


1.2+0.9 


PKS 0537-441 


70201116/17/18 


1997 Oct 18 


PKS 0537-441 


18.4±2.0 


29.1+5.5 


32.6+5.3 


0.8+0.2 


PKS 0735+178 


72301882/83/84 


1997 Nov 8 


PKS 0735+178 


19.3±2.0 


23.2+6.4 


33.4+5.1 


0.7+0.3 


PKS 1 144-379 


25400637/38/39 


1996 Jul 27 


PKS 1144-379 


10.7+1.7 


14.4+6.2 


37.2+5.7 


1.6+0.3 


OQ 530 


40200746/47/48 


1996 Dec 22 


OQ 530 


21.0+2.6 


32.8+5.3 


44.1+5.2 


1.0+0.2 


PKS 1519-273 


43401152/53/54 


1997 Jan 23 


PKS 1519-273 


9.9±2.9 


11.2+6.4 


<14.1 


0.4+1.6 


PKS 1749+701 


39901458/59/60 


1996 Dec 19 


PKS 1749+701 


29.2±2.5 


46.7+6.6 


54.3+5.5 


0.8+0.2 


PKS 1749+096 


82200287/88/89 


1998 Feb 14 


PKS 1749+096 


33.2±2.4 


43.0+7.3 


68.8+6.2 


1.0+0.2 


S5 1803+784 


15100861/62/63 


1996 Apr 16 


S5 1803+784 


29.3+2.0 


33.0+5.7 


66.8+5.1 


1.1+0.1 


4C 56.27 


15100964/65/66 


1996 Apr 16 


4C 56.27 


6.1+2.2 


8.8+3.7 


9.2+3.7 


0.5+0.6 


PKS 2131-021 


38601073/74/75 


1996 Dec 6 


PKS 2131-021 


12.1+2.3 


<17.6 


10.7+4.3 


-0.2+0.5 


PKS 2240-260 


16402576/77/78 


1996 Apr 29 


PKS 2240-260 


6.9+2.1 


21.6+5.9 


17.8+5.3 


1.2+0.5 


PKS 2254+074 


37900579/80/81 


1996 Nov 29 


PKS 2254+074 


14.1+2.3 


17.8+6.1 


30.5+5.7 


1.0+0.4 



"mJy 



spectively. A journal of the observations is presented in Tab.|2l 
which gives the name of the source, the corresponding ISOJd 
which uniquely identifies an ISO observation, and the observ- 
ing date. 

The data analysis was performed with the ISOCAM 
Interactive Analysis (CIA; Ott et al. ,1997j tool, Ver 5.0. The 
default steps of CIA were used for dark current subtraction, re- 
moval of cosmic ray hits (glitches), exposure co-addition, flat 
fielding, and flux conversion to astronomical units^. The tran- 
sient correction we adopted was not the standard one but that 
based on the Fouks-Schubert model (Fouks & Schubert '1995) 
which is valid for low contrasted ISOCAM LW observations 
(Coulais & Abergei r2000t . This was necessary since our LW 
observations were preceded by short-wavelength (SW) mea- 
surements (which could not be used), which made more diffi- 
cult to get the history of the LW detector. 

The source flux was derived by performing aperture pho- 
tometry with a radius which contained ~ 90% of the flux for 
each filter, namely 2 pixels for LW2, 3 pixels for LW7, and 4 
pixels for LW3, correcting for the missing flux using the CIA 
point spread function (PSF) library. The background was de- 
rived in a four pixel wide annulus whose inner radius was two 
pixels larger than the extraction radius for most sources (see 
below). The aperture centre was determined in an iterative way 
by first placing it on the brightest pixel and then by moving it 
in 1/4 pixel steps around that position until the largest value 
was reached. The flux uncertainty was estimated by adding in 
quadrature the root mean squared (rms) of the image within the 
extraction radius and the Icr uncertainty of the background es- 
timate. Systematic errors of ISOCAM photometry are ~ 15% 



This last step assumes oc v"', which is the average spectrum of 
our sources (see below), so no colour-correction was applied (in the 
most extreme case, fy oc v"--^, this would have been ~ 4%, way below 
our uncertainties). 



for point sources (Blommaert et al. l2003t . which are typically 
below our uncertainties. 

3.1 .1 . Notes on individual sources 

PKS 0138-097: faint source with a stronger source ~ 6.5 
pixels away in the LW2 image; we then estimated the back- 
ground in a two (instead of four) pixel wide annulus. 

PKS 0426-380: faint source; to increase the signal to noise 
(S/N) we used an extraction radius of 1 pixel for the LW2 filter, 
instead of our default value of two. No significant detection in 
the LW7 filter. 

S5 0454+844: faint source with another source of compara- 
ble flux ~ 8 pixels away in the LW2 image; we then estimated 
the background in a three (instead of four) pixel wide annulus. 
No significant detection in the LW7 filter Faint source with an- 
other source of comparable flux ~ 8 pixels away in the LW3 
image; we then estimated the background in a four pixel wide 
annulus just outside the extraction circle. 

PKS 0735+178: another source of comparable flux ~ 5.5 
pixels away in the LW7 image; we then estimated the back- 
ground in a four pixel wide annulus whose inner radius was 
four (instead of two) pixels larger than the extraction radius. 
Same for the LW3 filter (another source ~ 6.5 pixels away). 

PKS 1519-273: another source of comparable flux ~ 7.5 
pixels away in the LW7 image; we then estimated the back- 
ground in a two (instead of four) pixel wide annulus. No sig- 
nificant detection in the LW3 filter. 

PKS 1749+701: structure in the background of the LW3 
image; we then estimated the background in a region closer to 
the source (inner radius was one pixel larger than the extraction 
radius and annulus was three pixels wide). 

4C 56.27: faint source; to increase the signal to noise (S/N) 
we used an extraction radius of three pixels for the LW3 filter, 
instead of four. 
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PKS 2131-021: faint source; no significant detection in the 
LW7 filter. To increase the signal to noise (S/N) we used an 
extraction radius of three pixels for the LW3 filter, instead of 
four 

PKS 2240-260: structure in the background of the LW2 
image; we then estimated the background in a region closer to 
the source (inner radius was one pixel larger than the extraction 
radius). 

3.2. ISOPHOT 

All sources were observed with ISOPHOT (Laureijs et al. 
[2003 1 in rectangular chopped mode, i.e., the line-of-sight was 
switched periodically between the source and the background 
position, the latter being placed at 180" from the source. For 
each source, data were taken with the P2 detector at 25//m, with 
the ClOO detector at 65, 80 and lOOyum, and with the C200 de- 
tector in the 120, 170 and 180;um bands. Typical on-source ex- 
posure times were 128 s (P2), 32-64 s (ClOO), and 32 s (C200). 
A journal of the observations is presented in Tab.|2] which gives 
the name of the source, the corresponding ISOjd, and the ob- 
serving date. 

The data reduction was performed using the ISOPHOT 
Interactive Analysis Software Package VI 0.0 (PIA, Gabriel et 
al. [T998j . After corrections for non-linearities of the integration 
ramps, an 8-point signal pattern was created from each obser- 
vation by overplotting and averaging the basic blocks of the 
observation (the repeated background-nsource cycles; for de- 
tails see Abraham et al. l2003> . The signals of the patterns were 
transformed to a standard reset interval of 1/4 s, and an orbital 
dependent dark current was subtracted. The signals were cor- 
rected for non-linearities of the detectors by applying signal 
linearization corrections. From the patterns a (source - back- 
ground) difference signal was then extracted. In most cases the 
measured difference signal underestimates the real signal due 
to short term detector transients. Therefore, we applied the sig- 
nal loss correction algorithm as implemented in PIA. This cor- 
rection is a function of the chopping frequency and the mea- 
sured difference signal. The flux calibration of the C100/C200 
measurements was performed by comparison with the on-board 
fine calibration source (FCS), which was also measured in 
chopped mode. At 25 yum the detector's actual sensitivity could 
be reliably predicted from the orbital position of the observa- 
tion, and an orbital dependent default responsivity was applied. 
Finally, the derived flux densities were corrected for the finite 
size of the aperture by using the standard correction values. 

In a chopped observation the slow detector baseline varia- 
tions are expected to be cancelled via the frequent alternation 
of the source and background positions. In the case of very 
faint sources, however, the time per chopper plateau was often 
rather long, and in many cases only two ON-OFF cycles were 
performed. In these observations the baseline variations may 
not be perfectly corrected via the chopping procedure, since the 
temporal evolution of the baseline was only scarcely sampled. 
From a comprehensive calibration study we found that in the 
case of P2 observations of faint sources with only two chopper 
cycles the resulting fluxes were incorrect, usually overestimat- 



ing the true value. In order to cancel the signal transient, we 
introduced a new correction at the beginning of the process- 
ing scheme, by applying a transient correction developed for 
staring (rather than chopped) P2 observations. 

The staring transient correction normally works on signals 
at the Signal per Ramp Data (SRD) level, but in the case of 
a chopped observation it has to be applied before the pat- 
tern is created. Thus we modified the PIA procedure "pro- 
cess_erd_to_pattern.pro" and applied the transient correction on 
each pair-wise signal before pattern creation. The agreement 
with the fluxes predicted by extrapolating the ISOCAM fluxes 
was significantly improved with the new correction. We applied 
this algorithm to all our P2 observations which had only two 
ON-OFF cycles. We then processed further the measurements 
as normal chopped observations. Notice that some residual 
problems with uncorrected P2 data might remain (see § 14.2. H . 

The flux uncertainty at 25 /im was estimated to be of the 
order of 50mJy at faint flux levels. For the C100/C200 filters 
one important source of systematic uncertainty may be related 
to the signal loss correction. In addition, the actual responsiv- 
ity of the detector had to be derived from the accompanying 
FCS measurement, thereby introducing an additional uncer- 
tainty factor. Finally, at long wavelengths (A > 100/jm) sky 
confusion noise can be the dominant source of photometric er- 
ror. 

Note that for the four sources for which ISOCAM detected 
a close-by, relatively strong source, namely PKS 0138-097, S5 
0454-H844, PKS 0735+178, and PKS 1519- 273, so me confu- 
sion is expected in the ISOPHOT band (see § I4.2.H . 

3.3. Results 

Our ISOCAM fluxes are reported in Tab. |3l When no signif- 
icant detection was obtained we give 3cr upper limits. Tab. |3] 
gives also the best fit spectral slopes (in frequency space) to the 
ISOCAM data. These were derived by taking into account the 
errors on the fluxes and by excluding non detections. The mean 
value is (oy) = 0.95 + 0.15 (median 1.0), while the weighted 
mean is (ay) - 0.99+0.06. This is consistent with expectations, 
as a synchrotron peak frequency in the IR means vfy x const 
in the ISO bands, and therefore fy oc v as indeed observed. 

Our ISOPHOT fluxes are reported in Tab. 0] When no sig- 
nificant detection was obtained, or the formal uncertainty was 
larger than the measured flux density value, 3cr upper limits are 
given. 

4. Discussion 

4.1. Host galaxy contribution to the IR flux 

Before studying the spectral energy distributions (SEDs) of our 
sources it is important to assess a possible contribution from 
the galaxy which hosts the BL Lac. To this aim, we estimated 
the host galaxy flux in the IR band, following the approach 
of Bertone et al. (2000). A literature search actually revealed 
that 13 out of our 17 BL Lacs are optically unresolved, that 
is no host galaxy is detected (Pursimo et al. 120021 Sbarufatti 



Table 4. ISOPHOT fluxes^ 
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Name 


P2/25 
[25/im] 


C100/C.50 
[65/jm] 


C100/C_70 
[80yum] 


C100/C_100 
[lOOytim] 


C200/C_120 
[120yum] 


C200/C_160 
[170yum] 


C200/C.180 
[180yum] 


PKS 0048-097 


195±48 


218+115 


<552 


<378 


<1032 


<1377 


< 933 


PKS 0118-272 


53±48 


196+69 


290+174 


150+132 


< 948 


331+249 


< 822 


PKS 0138-097 


199±48 


317+76 


271+65 


<549 


<1089 


<1524 


< 981 


PKS 0426-380 


63±48 


126+60 


<819 


<180 


< 792 


400+201 


< 762 


S5 0454+844 


<144 


<375 


<783 


<582 


<1746 


<3771 


<6423 


PKS 0537-441 


70±48 


249+60 


155+117 


<180 


365+315 


404+387 


< 933 


DVQ n7QCi 170 


Zzj±4o 




<jOO 


1 A"? 
j23+lo/ 


^ 1 1 Q7 


^ Q7'^ 


<LVZ.y 


PKS 1 144-379 


115±48 


269+75 


259+61 


<309 


<1353 


<1566 


<1857 


OQ 530 


139±48 


160+97 


183+60 


204+181 


< 837 


351+237 


337+251 


PKS 1519-273 


249±48 


312+77 


164+60 


137+115 


<1305 


<4215 


<5025 


PKS 1749+701 


132±48 


<285 


76+60 


247+60 


< 846 


<2307 


<1725 


PKS 1749+096 


134±48 


503+75 


619+93 


346+102 


<1347 


<3495 


<4146 


S5 1803+784 


162±48 


75+60 


330+274 


312+60 


410+305 


472+278 


507+300 


4C 56.27 


63±48 


712+144 


<429 


<411 


<1182 


<1914 


459+416 


PKS 2131-021 


178±48 


<189 


<282 


170+60 


<1017 


< 927 


387+331 


PKS 2240-260 


139±48 


<195 


<213 


198+60 


< 909 


<1146 


<1122 


PKS 2254+074 


190+48 


114+96 


379+115 


265+185 


<1017 


<1104 


<1134 



"mJy 



et al. |2p5\. For the remaining four sources, we found host 
galaxy R magnitudes in Pursimo et al. (2002) and converted 
them to B magnitudes using the evolutionary synthesis model 
of Dunlop et al. ( 1989 1. We then used the results of Mazzei & 
DeZotti ( 1994 1, who calculated the flux ratio between the IRAS 
and the B bands for a sample of 47 ellipticals, to estimate the 
galaxy contributions at 12, 25, 60, and 100 fim. These turned 
out to be extremely low. The maximum expected contributions 
in fact, where around 1%, 0.01%, 0.1%, and 0.2% of the ob- 
served fluxes at 12, 25, 60, and 100 jim respectively. We thus 
conclude that for the BL Lacs in our sample the contribution of 
the host galaxy to their ISO flux is completely negligible. 

4.2. Spectral energy distributions 

To constrain the synchrotron power peak and to address the 
relevance of our ISO data in terms of emission processes in 
BL Lacs, we have assembled multifrequency data for all our 
sources. We have also included S5 2007+777, a 1 Jy BL Lac at 
z = 0.342 with previously published ISOPHOT data (Peng et 
al.|2000l. Since we are not interested in characterizing the in- 
verse Compton emission, we limited ourselves to frequencies < 
10'^ Hz. The main source of information was NED. Additional 
data were taken from the WMAP catalogue (Bennett et al. 
I2OOT1, the Two Micron All Sky Survey (2MASS) (Cutri et al. 
I2OOO1, the Guide Star Catalogue II'*, and the Sloan Digital Sky 
Survey (SDSS). Therefore, most data are not simultaneous with 
our observations. For more than half of our sources, however, 
we were able to find nearly-simultaneous (typically within a 
month) radio observations in the University of Michigan Radio 
Astronomy Observatory (UMRAO) database^. These are re- 

htpp: //archive. stsci.edu/gsc/ 
' The Bologna group has also performed a radio monitoring of 
some of our sources (Venturi et al. l2001> . The UMRAO observations, 
however, were always closer in time to our observing dates. 



ported in Table |5] Furthermore, nearly-simultaneous optical 
data were available f or PK S 0735+178 (Bai et aI. I999 l, OQ 
530 (Massaro et al. l2004l . and S5 1803+784 (Nesci et al. 
■^OOTl. As regards S5 2007+777, Peng et al. ('2000'l published 
simultaneous radio (seven frequencies), ISOPHOT (60 and 100 
fim), and optical (R band) observations of this object. 

The SEDs for our sources are shown in Fig.^ where big- 
ger, darker symbols and upper limits indicate ISO data and the 
nearly-simultaneous radio and optical data, and smaller, lighter 
symbols represent non-simultaneous data collected from the 
literature. 

4.2.1 . Synchrotron peak frequencies 

We determined Vpeak for our sources by applying an homo- 
geneous synchrotron - inverse self-Compton (SSC) model in 
the log V - log v/v plane to the SEDs. The SSC model was 
adapted from Tavecchio et al. ( 1998 1 and assumes that radia- 
tion is produced by a population of relativistic electrons emit- 
ting synchrotron radiation in a single zone of a jet that is mov- 
ing at relativistic speed and at a small angle to the line of sight 
(see Giommi et al. 2002 and Padovani et al. 120031 for previ- 
ous applications of this model and more details). Padovani et 
al. ( i2003i) make the point that, despite the fact that we can- 
not fully constrain the model parameters, Vpeak, a combination 
of magnetic field strength, Doppler factor, and electron break 
energy, is relatively well constrained even for relatively large 
variations of these parameters, which affect more strongly the 
high-energy/inverse Compton part of the SED. 

Some previous papers have derived Vpeak by fitting analyti- 
cal functions, such as a parabola or a third-degree polynomial, 
to the SED of blazars (e.g., Sambruna et al. 1996 Fossati et 
al. I1998> . We believe that our approach, although more com- 
plex and time consuming, is more robust especially when deal- 
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Table 5. Nearly-simultaneous radio observations. 



Name F4 sghz Observing date Fg ooHz Observing date Fi4.5ghz Observing date 

(Jy) (Jy) (Jy) 



PKS 0048-097 


1.58 = 


t0.06 


1996 


Dec 


13 


1.96 = 


t0.03 


1996 


Dec 9 


2.12 = 


t0.05 


1996 


Dec 11 


PKS 0735+178 


1.06 = 


t0.04 


1997 


Dec 


2 


0.69 = 


t0.16 


1997 


Nov 24 


0.94 = 


tO.ll 


1997 


Dec 10 


OQ 530 












0.56 = 


t0.15 


1997 


Jan 6 


0.64 = 


t0.03 


1996 


Dec 30 


PKS 1749+701 


0.63 = 


t0.04 


1996 


Dec 


13 










0.86 = 


t0.03 


1996 


Dec 19 


PKS 1749+096 


2.26 = 


t0.03 


1998 


Feb 


18 


3.40 = 


t0.08 


1998 


Feb 17 


4.02 = 


t0.09 


1998 


Feb 4 


S5 1803+784 


2.91 = 


t0.04 


1996 


Apr 


3 


2.90 = 


t0.19 


1996 


Apr 23 


3.18 = 


t0.06 


1996 


Apr 11 


4C 56.27 


1.79 = 


t0.04 


1996 


Apr 


16 


2.34 = 


t0.06 


1996 


Apr 24 










PKS 2131-021 












1.67 = 


t0.02 


1997 


Jan 6 


1.56 = 


t0.04 


1996 


Dec 19 


PKS 2254+074 




















0.48 = 


t0.02 


1996 


Dec 19 



ing with sparsely sampled SEDs, as we are guided by physics 
rather than just analytical fitting. 

During our fitting procedure we adhered to the follow- 
ing guidelines: 1. more weight was given to our nearly- 
simultaneous data; 2. more weight was given to ISOCAM data, 
which have smaller error bars; 3. clear outliers were not taken 
into account; in case of ISOPHOT data, these indicate confu- 
sion (see § l3.2t and/or some residual calibration issues; 4. only 
data points with v ^ 3 x 1 Hz were included. Our derived rest- 
frame Vpeak values are given in Tab.|6l Based on our experience, 
we beheve that most of our values are determined to better than 
0.5 dex. Exceptions to this rule include PKS 1519-273, PKS 
2131-021, and PKS 2254+074, whose SEDs are sparse and/or 
whose ISO continuum shape is somewhat iiTeguIar (see Fig.^. 

The Vpeak distribution is very narrow, with 10/18 sources 
having values in the range 1 - 3 x 10'^ Hz. The mean value is 
(logVpeak) - 13.0 + 0.1 (median 13.0), that is ~ 30yum. Given 
our set of simultaneous IR data, these represent the best deter- 
minations available of the synchrotron peak frequencies for the 
1 Jy BL Lac sample. 



4.2.2. Synchrotron peak variability in IR band 

Inspection of Fig. \l\ shows that quite a few of our sources 
have a nearly-simultaneous SED around the peak which is 
quite different from the one previously available using non- 
simultaneous data from the literature. The most obvious cases 
include PKS 0048-097, PKS 0537-441, PKS 0735+178, OQ 
530, 4C 56.27, where Vpeak could have changed by over one 
order of magnitude. This impression is confirmed by a com- 
parison with the Vpeak values derived by Sambruna et al. ( 1996 1 
for these objects, which are larger than ours by factors ~ 20, 
120, 10, 2, and 8 respectively. This may be indicative of large 
changes in the location of the synchrotron peak, as observed 
in the case of the high-energy peaked BL Lac MKN 501 (e.g.. 
Plan et al. 199SL Massaro et al. 2006 1. Alternatively, Vpeak could 
have stayed constant, and the large optical variability would 
be due to a shift in normalization of the whole SED, imply- 
ing a coiTespondingly high change in the IR flux. Simultaneous 
optical-IR monitoring is required to be able to distinguish be- 
tween the two possibilities. 



Table 6. Rest-frame synchrotron peak frequencies and wave- 
lenghts. 



Name 


Vpeak 






10'3 Hz 


yum 


PKS 0048-097 


0.3 


100 


PKS 0118-272 


10 


3 


PKS 0138-097 


0.3 


100 


PKS 0426-380 


0.5 


60 


S5 0454+844 


0.5 


60 


PKS 0537-441 


0.1 


300 


PKS 0735+178 


1 


30 


PKS 1 144-379 


1 


30 


OQ530 


3 


10 


PKS 1519-273 


1 


30 


PKS 1749+701 


3 


10 


PKS 1749+096 


0.3 


100 


S5 1803+784 


1 


30 


4C 56.27 


0.5 


60 


S5 2007+777" 


1 


30 


PKS 2131-021 


2 


15 


PKS 2240-260 


2 


15 


PKS 2254+074 


2 


15 


"ISOPHOT data published by Peng et al. 


l2000t. 



4.3. Astrophysical relevance 

A precise measurement of Vpeak can obviously characterize syn- 
chrotron emission from blazars. Furthermore, its value and that 
of the spectral curvature around the maximum emission can 
constrain particle acceleration mechanisms in these sources 
(e.g., Massaro et al. 2006 1. This is also extremely relevant for 
two up-coming y-ray missions, that is AGILE and GLAST 
These observatories, in fact, will sample inverse Compton 
emission above the peak for LBL and flat-spectrum radio 
quasars. Moreover, the large synchrotron peak variability pos- 
sibly present in a few of our sources will also be important to 
constrain the link between the two main emission processes in 
blazars and the source of the inverse Compton seed photons 
(e.g., Sokolov et al. .2004j . 

4.4. Previous ISO observations of BL Lacs 

To the best of our knowledge, only three ISO observations of 
BL Lacs have been published. Beside the already mentioned 
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S5 2007+777, these include PKS 2155-304 (ISOCAM and 
ISOPHOT; Bertone et al. 2000 1 and Mrk 180 (ISOPHOT, one 
detection at 90yum; Anton et al. i2004j . both of the HBL type. 

5. Conclusions 

We have presented new, simultaneous ISOCAM and ISOPHOT 
observations covering the 7-200/im range for seventeen (out of 
34) BL Lacertae objects belonging to the 1 Jy sample, the pro- 
totype radio-selected (low-energy peaked) sample. These data 
open up a new window on the spectral energy distributions of 
low-energy peaked BL Lacs and fill a glaring gap in their far- 
infrared spectra. They also represent a sevenfold (eighteenfold) 
increase in the number of published ISO results for BL Lacs 
(LBL). 

Our main results can be summarized as follows: 

1. Our detection rate is good; 100%, 82%, and 94% at 6.7 
(LW2), 9.6 (LW7), and 14.3 (LW3) //m (ISOCAM), and 
94%, ~ 65%, and ~ 20% at 25 (P2), 65 - 100 (ClOO), 
and 120-180 (C200) fim (ISOPHOT) respectively. To deal 
with the effect of chopping on our relatively faint sources, 
which would lead to a flux overestimate at 25 ;um, we have 
also employed a novel correction in the ISOPHOT data re- 
duction; 

2. We have built the spectral energy distributions of our 
sources, with the addition of a previously published one, 
by using our new ISO data, complemented by nearly- 
simultaneous radio data for about half the sample (and op- 
tical data for four sources), and other non-simultaneous 
multi-frequency data. We have then derived the synchrotron 
peak frequency, Vpeak- This turned out to have a very nar- 
row distribution, with ~ 60% of the sources having val- 
ues in the range 1 - 3 x 10'^ Hz, centered at ~ 10'^ Hz 
(~ 30;um). These represent the best determinations avail- 
able of the synchrotron peak frequencies for low-energy 
peaked BL Lacs; 

3. A comparison with previous such estimates, based on non- 
simultaneous optical and near infrared data, may indicate 
strong Vpeak variations in a number of sources, possibly as- 
sociated with large flares as observed in the high-energy 
peaked BL Lac MKN 501. 

These results, apart from characterizing synchrotron emis- 
sion in blazars, can constrain particle acceleration mechanisms 
and are also important for the planning and data interpretation 
of up-coming new y-ray missions. 
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Fig. 1. The spectral energy distributions for our sources. Bigger, darker symbols indicate ISO data and the nearly-simultaneous 
radio (and optical in the case of PKS 0735+178, OQ 530, S5 1803+784, and S5 2007+777) data (arrows denote ISO upper limits), 
while smaller, lighter symbols represent non-simultaneous data collected from the hterature. The dashed line is our synchrotron 
model. 
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Fig. 1. (continued) The spectral energy distributions for our sources. Bigger, darker symbols indicate ISO data and the nearly- 
simultaneous radio (and optical in the case of PKS 0735+178, OQ 530, S5 1803+784, and S5 2007+777) data (arrows denote 
ISO upper limits), while smaller, lighter symbols represent non-simultaneous data collected from the literature. The dashed line 
is our synchrotron model. 
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Fig. 1. (continued) The spectral energy distributions for our sources. Bigger, darker symbols indicate ISO data and the nearly- 
simultaneous radio (and optical in the case of PKS 0735+178, OQ 530, S5 1803+784, and S5 2007+777) data (arrows denote 
ISO upper limits), while smaller, lighter symbols represent non-simultaneous data collected from the literature. The dashed line 
is our synchrotron model. 
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Fig. 1. (continued) The spectral energy distributions for our sources. Bigger, darker symbols indicate ISO data and the nearly- 
simultaneous radio (and optical in the case of PKS 0735+178, OQ 530, S5 1803+784, and S5 2007+777) data (arrows denote 
ISO upper limits), while smaller, lighter symbols represent non-simultaneous data collected from the literature. The dashed line 
is our synchrotron model. 
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Fig. 1. (continued) The spectral energy distributions for our sources. Bigger, darker symbols indicate ISO data and the nearly- 
simultaneous radio (and optical in the case of PKS 0735+178, OQ 530, S5 1803+784, and S5 2007+777) data (arrows denote 
ISO upper limits), while smaller, lighter symbols represent non-simultaneous data collected from the literature. The dashed line 
is our synchrotron model. 



